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Context. In this paper we present the first SiO multiline analysis (from 7=2-1 to 7=11-10) of the molecular bullets along the outflows of the 
Class sources L1448-mm and L1157-mm, obtained through observations with IRAM and JCMT. 

Aims. This analysis has been performed to investigate the physical properties of different types of bullets, in order to explore their excitation 
conditions and put constraints on their formation and evolution. 

Methods. For this purpose we have computed the main physical parameters {n H2 , T kin , N S io) in each bullet and compared them with other 
tracers of warm and dense gas and with models for the SiO excitation in shocks. 

Results. We find that the bullets close to L1448-mm, associated with high velocity gas, have higher excitation conditions {n Hl ~10 6 cm -3 , 
T k, 500 K) with respect to the LI 157 bullets ( n Hl ~l-5xl0 5 cm -3 , T ~ 100 - 300 K). In both the sources, there is a clear evidence of the 
presence of velocity components having different excitation conditions, with the denser and/or warmer gas associated with the gas at the higher 
speed. In L1448 the bulk of the emission is due to the high-excitation and high velocity gas, while in LI 157 most of the emission comes from 
the low excitation gas at ambient velocity. The observed velocity- averaged line ratios are well reproduced by shocks with speeds v s larger than 
~30 km s _1 and densities ~ 10 5 - 10 6 cm -3 . Plane-parallel shock models, however, fail to predict all the observed line profiles and in particular 
the very similar profiles shown by both low and high excitation lines. 

Conclusions. The overall observations support the idea that the LI 157 clumps are shock interaction events older than the LI 448 bullets close to 
the driving source. In the latter objects, the velocity structure and the variations of physical parameters with the velocity resemble very closely 
those found in optical/IR jets near the protostar, suggesting that similar launching and excitation mechanisms are also at the origin of collimated 
jets seen at millimetre wavelengths. 
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ABSTRACT 



1. Introduction 



of an embedded atomic jet, or whether the jet itself is intrin- 
sically molecular. In particular, the EHV molecular jets show 
a knotty structure with regularly spaced clumps, the so-called 
molecular bullets, that are well-defined entities travelling at 
velocities of hundreds of km s" 1 (e.g. Bachiller et al. 119901 
Chandler & Richer 2001 ). The bullets are most likely associ- 
ated with shocks formed by the outflow propagation into a ho- 
mogeneous material and are consequently thought to be a di- 
rect manifestations of the episodic nature of outflows (Dutrey 
et al. fTW7l Arce & Goodman I200T1 120021) . On the other hand, 
a few outflows driven by Class objects (e.g. L1157-mm, 
NGC2264G, BHR71) show another type of molecular clump, 
characterised by lower velocities and by a rich chemistry: these 
clumps have strong emission due to molecular species which 
increase their abundance only when a shock-induced chemistry 
is active (Bourke et al. 119971 Bachiller et al. 2001 ). The origin 
of these chemically rich clumps is not clear: a possibility is that 



The mass-loss process starts during the earliest phases of star 
formation when protostars are still accreting mass (e.g. Richer 




that interact with the high-density clump hosting the newly- 
born stars, forming shocks and driving molecular outflows. 
Fast (> 100 km s" 1 ) collimated jets are observed both in op- 
tical emission of atomic lines and near infrared (NIR) H2 emis- 
sion, showing excitation conditions with T ex ranging from ^ 
2000 K to more than 10000 K. Also, extremely high veloc- 
ity (EHV) jets have been seen in the form of narrow and 
high velocity CO and SiO molecular outflows in a handful of 
Class sources (e.g. L1448-mm, HH211, IRAS04166, Gueth 
& Guilloteau 1999, Tafalla et al. 2004). It is a matter of some 
debate whether the EHV molecular jet is the cold external layer 
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they are a consequence of a pre-existing density enhancement 
and of its interaction with the passage of a shock (Viti et al. 

However, the nature of both kinds of clumps and the possi- 
ble relationships between them based on age, physical charac- 
teristics, and/or excitation conditions, are not yet understood. 
Detailed quantitative analysis of the EHV jets, bullets, and 
chemically rich clumps have been so far obtained with two 
standard outflow tracers such as CO and SiO, but only with low 
frequency and consequently low excitation transitions. These 
transitions, however, are not suited to probe the warm gas com- 
ponents that may be associated with high- velocity and shocked 
gas. Indeed, sub-mm- wavelengths observations show that some 
molecular bullets are dense, n Hl ~ 10 6 cm -3 , and warm, T ex 
> 100 K (Hatchell et al. IT9971 Nisini et al. I5UU21 Palau et al. 
120061 

In addition, ISO-LWS high-J CO and H2O low angular res- 
olution observations seem to suggest that EHV jets are associ- 
ated with temperatures in the 300-1500 K range (Giannini et 
al. 120011) . Whether this warm gas represents a major compo- 
nent of the mass flux and consequently of the energetics of the 
mass loss from a YSO is still an open question. 

In this paper, we investigate the nature of both the high- 
velocity bullets and chemically-rich clumps through the use of 
a multi-frequency analysis of the SiO emission. This molecule 
is the most typical tracer of high-density conditions and high 
velocity shocks, strong enough to release refractory elements 
in the gas phase (e.g. van Dishoeck & Blake 1998 ). So far, SiO 
has been mainly observed at high velocities through its low 
rotational transitions (/ < 5), tracing relatively cold (< 100 K) 
gas (M artin-Pintado et al. [T^92l Codella et al. [T999l Gibb et al. 
120041 

The objects studied in this paper are the outflows LI 44 8 
and LI 157, representative of flows associated with EHV jet- 
like bullets and lower velocity clumps respectively. In partic- 
ular, here we complement the low excitation spectra of SiO 
7=2-1, 3-2 and 5-4, with new observations of the / = 6-5, 
8-7, and 11-10 emission at submillimetre wavelengths. High- 
/ SiO emission from / = 5-4 to 11-10 has been already ob- 
served in another typical EHV flow, HH211, by Nisini et al. 
(2002). Here we perform a detailed analysis of the full set of 
lines with the following aims: (i) disentangle, through their ex- 
citation conditions, the warm gas components directly origi- 
nating in the EHV jet from ambient and colder gas entrained 
by the jet itself; (ii) compare the spectra due to different tran- 
sitions and thus investigate whether the excitation conditions 
vary with velocity, comparing the results with the predictions 
given by shock models; and (iii) estimate the type and strength 
of the shocks, and so place constraints on the origin of bullets 
and chemically rich clumps. 

2. Observations and results 

Observations were carried out with the IRAM 30-m telescope 
for the SiO transitions from / = 2-1 to 6-5, and with the James 
Clerk Maxwell Telescope (JCMT) for the SiO transitions from 
/ = 8-7 to 11-10. In Table 1, a list of the lines observed and 



their excitation energy is presented together with relevant ob- 
servational parameters. 

JCMT observations were performed during several nights 
in November 2001, and October and November 2002. We used 
the heterodyne receiver RxW for the / = 11-10 and / = 10-9 
observations (HPBW = 12"), and RxB3 for the / = 8 - 7 
observations at 347 GHz (HPBW = 14"). The back-end was 
the digital autocorrelation spectrometer (DAS) with a band- 
width of 500 MHz. To convert the antenna temperature (7^) 
into main-beam brightness temperatures, main-beam efficien- 
cies t]mb equal to 0.55 and 0.4 were adopted for receivers B3 
and W respectively. Pointing was checked on NGC7538 IRS1 
and W3(OH) and was good within about 3". 

In LI 448, we observed strips of 5 and 3 points in the SiO 
8-7 and 11-10 lines respectively, along the major axis of the 
outflow (PA. 157?4) with a 10" spacing (see FigQ}. The (0,0) 
position of this strip corresponds to the exciting source of the 
outflow, L1448-mm feooo=03 /z 22 m 34.3', £ 2 ooo=30 J 33 m 350. 

The pointings along the strip covers the CO bullets Rl 
and R2 identified by Bachiller et al. (1990), resolved by the 
Guilloteau et al. (1992) interferometric observations in two dis- 
tinct peaks each, named BI-BII and RI-RII and indicated in 
FigE 

In the (0,0) position a spectrum of the SiO 10-9 transition 
was also obtained. Spectra of the pointings centered on the 
mm source and on the B1,R1 bullets are shown in Fig. |4] 
In addition to the central region, a SiO 8-7 spectrum was 
obtained also in the south lobe of the outflow, at the position 
were a strong SiO 2-1 peak was observed by Dutrey et al. 
1997 (position R4 following Bachiller et al. (1990)). 

In LI 157, a map in the SiO 8-7 line was acquired covering 
the strongest bullets identified by Bachiller et al 2001 , B0, 
Bl and B2 in the southern lobe of the outflow, and R0, Rl 
and R in the northern lobe (Fig. 13 . Additional spectra of the 
SiO 10-9 transitions were obtained only in the B0, Bl and B2 
positions. Fig.[3]shows the spectra of the brightest Bl peak. 

The JCMT observations were complemented by observa- 
tions of SiO transitions at lower frequency obtained at nearby 
positions using the IRAM-30m telescope (Santiago Garcia, 
PhD thesis). These observations were performed on the 14th 
June 2000. Heterodyne receivers A100, A230, D150, and 
D270 were tuned to observe all the lines simultaneously in 
wobbler- switching mode, with HPBW beams from 28"to 
9", and a wobbler throw of ~240". For all the transitions an 
integration time of 4 minutes was accumulated in all the LI 157 
bullet spectra, and in the Bl, Rl and R4 spectra in LI 44 8, 
while the spectrum of the central L1448-mm position was inte- 
grated for 8 minutes. Pointing was checked in nearby sources 
and found to be accurate within 2". System temperatures 
ranged from ^166 K - 208 K for the J=2-l line, to « 1016 K 
- 2180 K for the J=6-5 line. Spectra obtained with IRAM are 
also shown in Figs. [3] and 01 together with the JCMT spectra. 

For both JCMT and IRAM observations, we have obtained 
the intensity integrated over the velocity channels, separating 
the blue- and red- shifted components. These values are listed in 
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Fig. 1. The pointings adopted for our observations of the cen- 
tral region of LI 448 are indicated by black squares over an 
interferometric velocity channel map at 64.8 km s" 1 of the SiO 
2-1 transition taken from Guilloteau et al. (1997). Offsets (in 
arcsec) are with respect to L1448-mm. 



Tables 2 and 3 only for the positions which will be relevant for 
our analysis. The errors on the integrated intensities reported 
in these tables are the statistical errors due to thermal noise 
on the individual spectra. Calibration accuracy is about 10% 
for lines with / <8, and 20% for the /=10-9 and 11-10 lines, 
and these uncertainties must be considered when estimating 
the total errors on observed quantities. In the LI 157 flow, the 
coordinates of the beam centers in the JCMT observations do 
not always coincide with those of the IRAM observations, in 
which case we report in Table 3 the intensities obtained from 
the spectra closer in space. The angular difference of these 
spectra is typically of the order of 2-3 ": in Fig.|2| the different 
displacement of the IRAM and JCMT pointings is indicated. 

Comparing the LI 157 SiO /= 8-7 map with the SiO /= 
3-2 and 5-4 maps in Bachiller et al. (2001), we see a similar 
brightness distribution: the strongest SiO emission is associated 
with the blue lobe, and in particular with bullet Bl, the one 
most chemically active along the flow. In our map, we resolve 
the peak of the B0 bullet that remained unresolved in the 5-4 
IRAM maps. We do not detect any SiO 8-7 emission in the 
Rl and R0 bullets. We also find a general agreement between 
the high-J SiO distribution and the distribution of the warm 
gas mapped through the pure rotational H2 0-0 transitions by 
ISOCAM (Cabrit et al. 1999). 

In LI 44 8 the situation is somehow different. The strongest 
high-J lines are detected close to the millimetre source, in the 
Bl and Rl bullets, while the R4 bullet in the southern lobe, 
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Fig. 2. Map of the SiO / = 8-7 integrated emission in the 
LI 157 outflow. Offsets are with respect to L1157-mm. First 
contour and steps are 0.75 K km s" 1 and 0.60 K km s" 1 for the 
blue and red bullets, respectively (which correspond to 3cr in 
each case). The circle represents the JCMT HPBW, while the 
small crosses indicate the observed positions. Filled squares in- 
dicate the pointings of the IRAM data considered in our analy- 
sis. 



although very strong in SiO 2-1, is much weaker in the high-/ 
lines. This behaviour was already pointed out by Bachiller et al. 
(1991) but this result is much clearer in our 8-7, 10-9 and 11- 
10 observations. In L1448 there is no H2 0-0 map to compare 
the morphology of the SiO emission with the occurrence of 
warm gas; however, we point out that the Rl and B 1 bullets are 
those associated with the region of the strongest high-J CO and 
H 2 emission observed by ISO-LWS (Nisini et al. 2000). 

The SiO 2-1 interferometric observations of Guilloteau et 
al. (1992) have shown that there is a velocity pattern along the 
SiO jet, with the velocity linearly increasing from BI-RI (at 3- 
4") to BI-RII (at -10") . We have evidence of this behaviour 
in the shapes of the lines at different frequency, which are ob- 
served with different beam sizes. In the observations centered 
on the mm source, the blue and red-shifted emission peaks are 
located at velocities which decrease going from the 2-1 (Vlsr 
= -55 and +63km s" 1 for the B 1 and Rl bullets respectively) to 
the 6-5 lines (Vlsr = -43.2 and +68.1kms _1 ), this latter hav- 
ing a beam size of only 9" and thus including only the slower 
BI-RI components. 
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Table 1. Observed SiO transitions and instrumental setting 



Line 




Instr. 


Freq. 


Res. 


HPBW 


IB 




K 




GHz 


km s _1 


// 




SiO J=2-l 


6.2 


IRAM 


86.85 


3.45 


27.6 


0.77 


SiO 7=3-2 


12.5 


IRAM 


130.27 


2.30 


18.4 


0.64 


SiO 7=5-4 


31.2 


IRAM 


217.10 


1.38 


11.0 


0.56 


SiO 7=6-5 


43.7 


IRAM 


260.52 


1.15 


9.2 


0.51 


SiO 7=8-7 


75.0 


JCMT 


347.33 


1.10 


14 


0.55 


SiO 7=10-9 


114.5 


JCMT 


434.12 


1.73 


11 


0.4 


SiO 7=11-10 


137.4 


JCMT 


477.50 


0.79 


11 


0.4 



Table 2. Observed SiO intensities in L1448-mm 



Offset* 


Bullet 


7=2-1 


7=3-2 


7=5-4 


7=6-5 
f T MB dV b 


7=8-7 


7=10-9 


7=11-10 


(-3.1, +9.5) 


Bl 


6.3+0.3 


14.0+0.4 


15.8+0.7 


16.7+1.2 


5.3+0.5 




<1.2 




Rl 


1.9+0.3 


3.2+0.4 


3.8+0.7 








<1.2 


(-2.3, +3.4) 


Bl 


7.0+0.4 


13.0+0.4 


15.8+0.7 


19.0+1.2 










Rl 


3.7+0.4 


7.7+0.4 


11.2+0.7 


12.7+1.2 








(0,0) 


Bl 


6.6+0.2 


12.5+0.4 


13.0+0.4 


14.6+0.6 


12.0+0.5 


6.7+0.9 


<1.2 




Rl 


3.8+0.2 


9.4+0.4 


16.5+0.4 


18.8+0.6 


11.6+0.5 


8.7+0.9 


4.5+0.9 


(+2.3, -3.4) 


Bl 


5.1+0.4 


6.8+0.4 


6.4+0.7 


3.5+1.1 










Rl 


4.4+0.4 


9.4+0.4 


13.5+0.7 


20.6+1.1 








(+3.1,-9.5) 


Bl 


2.6+0.4 


2.6+0.5 




1.4+1.2 


4.0+0.4 




<1.2 




Rl 


4.0+0.4 


10.0+0.5 


11.1+0.7 


11.5+1.2 


13.3+0.4 




6.7+0.9 


(+26,-128) 


R4 


13.1+0.3 


17.1+0.5 


10.7+0.5 


8.3+1.2 


3.1+0.4 







a Offset in arcsec with respect to L1448-mm (of 2 ooo=03 /7 22 m 34.3*, £ 20 oo=+30 J 33 m 35 5 ) 

b Flux errors refer to the rms noise in the individual spetra only. Total calibration uncertainty are of the order of 10% for the 7 <8 lines and 
20% for the 7=10 and 11 lines. 



3. Analysis of the SiO excitation 

For a quantitative analysis of the SiO excitation employing ra- 
tios among different lines, we need to correct the derived in- 
tensities for the different beams adopted for the observations. 
In LI 44 8 this correction can be significant, especially for the 
Bl and Rl knots that are located at 20" from the mm source. 
Each of these bullets is about 10" in length, but is not resolved 
in the transverse direction at a resolution of 2" (Guilloteau et 
al. 1992). Since we base our excitation analysis on the ratios 
between different lines, we have convolved all the observations 
along the strip to the beam of the 7=2-1 transition, in order to 
obtain corrected ratios on the Bl and Rl bullets. In LI 157 the 
problem seems to be less severe, since SiO IRAM maps of dif- 
ferent transitions show that the bullets emission region is large 
even for th e SiO 2-1 line taken with a 27" beam (Bachiller at 
al. 2001). In order to properly account for beam-filling effects, 
we have convolved the 8-7 data to same resolution as the 2-1 
beam. We find that the peak temperature of the line decreases 
by a factor from 1.3 in knot B 2 to 2.2 in knot R. These changes 
in the peak temperature are consistent with a size of about 18- 
22 arcsec in knots R, B0, and Bl and around 30 arcsec for knot 
B2. In the analysis we have considered both the case of uncor- 
rected line ratios and line ratios corrected assuming the above 
estimated sizes in order to evaluate the effect of beam dilution 
on the final results. 



3.1. Large Velocity Gradient analysis 

For the excitation analysis, we use an LVG model that consid- 
ers the first 20 SiO levels, with collision rates from Turner et 
al. (1992) up to 300 K, plus the extrapolations given by Schoier 
et al. (2005) for larger temperatures (up to 2000 K). In Fig. [3 
we plot the expected line brightness as a function of the rota- 
tional quantum number, for different values of the gas kinetic 
temperature and H2 density. While lines with J up < 5 do not 
strongly depend on the temperature but are sensitive to the 
density, high temperatures regimes (T>100 K) can be distin- 
guished from low temperature gas only through high- J (7 > 5) 
lines. The plots also show that the lines become insensitive to 
the temperature for T^500 K. These plots have been obtained 
using a ratio between column density and velocity dispersion 
(Nsjo/AV) of 10 13 cm" 2 km" 1 s. In this specific case, the maxi- 
mum line opacity is ~ 0.2. Opacity effects start to become sig- 
nificant, affecting the line ratios up to the point at which the 
lines become insensitive to the SiO abundance, for Ns/o/AV ~ 
5xl0 13 cm" 2 km" 1 s, i.e. for column densities of the order of 
5xl0 14 cm" 2 in gas with total velocity about lOkms" 1 . Since 
this value is significantly larger than the typical SiO column 
densities derived in shocked regions (i.e. Codella et al. 1999), 
we have assumed in the first instance that the lines are optically 
thin in order to derive temperature and density from the differ- 
ent line ratios. Large values of N^o/AV have been introduced 
only when a reasonable fit could not be obtained assuming 
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Table 3. Observed SiO intensities in LI 157 



Offset* 


bullet 


7=2-1 


7=3-2 


7=5-4 
f Tmb 


7=6-5 

dV b 


7=8-7 


7=10-9 


(-30, +100) 


Rl 


4.7+0.2 


4.7+0.4 


<0.5 


<0.7 






(-30, +140) 


R 


10.1+0.2 


13.4+0.4 


7.8+0.9 


2.0+1.2 






(-28.5, +137.5) 












1.7+0.3 




(+22.5, -64.5) 


Bl 


19.3+0.3 


28.9+0.4 


24.9+0.9 


18.8+1.2 






(+21.1,-62.5) 












8.7+0.4 


2.6+0.6 


(+21.1,-72.5) 












4.6+0.4 




(+22.5, -36) 


B0 


7.4+0.3 


9.5+0.4 


6.2+0.7 


4.4+1.2 






(+21.1,-42.5) 












2.5+0.3 


<0.8 


(+35,-100) 


B2 


12.7+0.3 


14.4+0.4 


6.5+0.7 


4.1 + 1.2 






(+31.1,-102.5) 












1.5+0.3 




(+41.1, -102.5) 












1.8+0.3 


< 1 


(+41.5,-110) 


B2 


12.5+0.3 


15.2+0.4 


8.6+0.7 


5.0+1.2 






(+41.1, -112.5) 












1.7+0.3 





a Offset in arcsec with respect to L1157-mm, (of 2 ooo=20 /l 39 m 06.2 5 , £ 2 ooo=68 J 02 m 15.9 5 ) 

b Flux errors refer to the rms noise in the individual spetra only. Total calibration uncertainty are of the order of 10% for the 7 <8 lines and 
20% for the 7=10 and 11 lines. 



optically-thin emission. Once T and nn 2 have been estimated, 
the SiO column density has been derived from the absolute in- 
tensity of the 7=5-4 line. We selected this line because is the 
one with the smaller beam having good S/N ratios. From the 
derived column density and the velocity dispersion measured 
on our spectra, we have checked a posteriori the validity of our 
assumption of low opacity. 

We have applied the LVG code on the intensities of each 
of the observed bullet, integrated in all the velocity channels, 
separating only the blue- shifted and red- shifted components. 
The results of the LVG analysis are discussed below for LI 448 
and LI 157. 

3.1.1. L1448 

Fig. |6| shows the observed integrated intensities, normalised 
to the 5-4 transitions, of the bullets sampled by our observa- 
tions, compared with the LVG model computations. Errors on 
line ratios take into account both the rms in the spectra and 
the calibration uncertainties. Assuming a velocity dispersion of 
20km s" 1 , i.e. comparable to the FWHM of the observed lines, 
the modelled line ratios are insensitive to the column density 
up to values 10 15 cm~ 2 . Since our observations do not fully 
sample the region covered by the larger 2-1 beam, the applied 
convolution may not properly takes into account all the emis- 
sion region. In addition to this, the lack of a good fit may sig- 
nal the presence of gradients in the physical parameters, with 
components at lower density/temperature encompassed by the 
large beams. Consequently, we adopted the best model that fits 
only the lines with 7 >5, which have been obtained through a 
more similar beam size and thus trace similar gas components. 
A density of uh 2 = 10 6 cm 3 is well constrained by the observed 
ratios, while any temperature T ^500 K gives an equally good 
fit. We have also considered models assuming a higher density 
and lower temperature to explore whether other pairs of values 



can fit the data equally well. We find that within the errors, the 
line ratios of knots Bl and Rl can be consistent with a temper- 
ature of 300 K (see Fig. 5). 

In the bullet R4 the best fit is obtained for T=200 K and 
?x# 2 =2.5xl0 5 cm 3 consistent with the lower excitation nature 
of this clump with respect to Bl and Rl. 

SiO column densities, determined from the absolute ob- 
served intensity of the 5-4 line, are of the order of 3-3.5xl0 13 
cm -2 in all the three bullets. Values derived for the Bl and 
Rl clumps, however, are probably just lower limits, since the 
emission region does not fill the entire 1 1 arcsec beam of the 
5-4 observation. If the emission region is confined to a jet-like 
knot about 2" in diameter (see e.g. Guilloteau 1992 ), the corre- 
sponding column density would be enhanced by about a factor 
of 4, thus being of the order of 10 14 cm -2 . 

3.1.2. L1157 

Fig. shows the observed line intensities, normalized to the 
7=5-4 transition, in the bullets R, Bl and B2, where we per- 
formed the LVG analysis. In the figure, filled circles represent 
line ratios not corrected for beam filling effects; open circles 
represent ratios corrected for beam filling, assuming the emis- 
sion sizes estimated by convolving the 7=8-7 observations to 
the 7=2-1 beam (see previous section). This correction signif- 
icantly affects only the 2-1 and 3-1 transitions, so we firstly 
fitted only the lines with 7 >5. In Bl and B2 the line ratios are 
consistent with a temperature of ~ 300 K, a density ~3xl0 5 
cm 3 and optically-thin lines. The bullet R, on the other hand, 
is consistent with colder gas at T -100 K. These fits reproduce 
rather well the 2-1 and 3-2 uncorrected ratios. We have also 
tried to reproduce all the beam-filling-corrected data points, ob- 
taining the fit indicated by dashed lines in Fig.0 Good fits are 
obtained also in this case, by adopting slighter lower tempera- 
tures of -50-200 K. Column densities are 2, 5 and 8xl0 13 cm~ 2 
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Fig. 4. Spectra of the SiO line emission taken in positions at the base of the L1448 outflow. Offsets (in arsec) of the Bl and Rl 
positions are with respect to the L1448-mm coordinates. 



at bullets R, B2 and Bl respectively. Therefore, the LI 157 
clumps appear less excited than the L1448 Bl and Rl bullets, 
confirming a dependence of the excitation conditions with the 
velocity of the bullets. 

3.2. Excitation as a function of velocity 

The above LVG analysis has been performed on intensities in- 
tegrated over all the observed velocity channels, thus giving 
physical parameters averaged over all the emitting gas. The 
observed line profiles, however, are not symmetric and clearly 



show that gas components at different velocities are seen along 
the line of sight. We can then investigate if these components 
are associated with different excitation conditions. 

Fig-Elshows that in the LI 157-B1 spectra there are at least 
two components, one peaking close to ambient velocity and 
one at about -15kms _1 , whose relative T B ratio increases go- 
ing from the 2-1 to the 8-7 line (see also Fig. 8). The signal- 
to-noise ratio in the line wings of the high-J lines is not high 
enough to perform a separate LVG analysis in the different ve- 
locity channels. We can however still infer a quantitative trend 
in the excitation as a function of velocity by comparing sensi- 
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Fig. 3. Observed SiO spectra towards the LI 157 Bl position. 

tive line ratios. In Fig. [8] we have plotted the SiO 8-7/5-4 ratio 
as a function of velocity for the LI 157 Bl bullet and for the 
L1448 Rl and Bl bullets. The 8-7 and 5-4 lines have been 
observed with similar beamsizes and should thus be not much 
affected by different beam fillings. 

For L1157-B1 the 8-7/5-4 ratio is -0.25+0.03 in the peak 
close to ambient Vlsr velocity and -0.6+0.25 at velocities 
15km s" 1 . In L1448 Bl/Rl, line ratio differences with ve- 
locity are less clear, although a trend is still recognisable in the 
Bl bullet, where the 8-7/5-4 ratio varies from -0.2-0.3+0.2 at 
Vlsr ~ -20km s" 1 and -0.8+0.08 at V LS r ~ -50km s" 1 . A ra- 
tio -0.2-0.3 is consistent with gas at a density of -10 6 cm 3 and 
temperature -200 K, while a value of 0.5-0.7 would require ei- 
ther a larger temperature (> 500 K) or a larger density (~5xl0 6 
cm 3 ). 

In L1448 Bl/Rl the observed excitation as a function of 
velocity is consistent with gas coming from different layers be- 
hind a shock: the maximum temperature is found at the shock 
front, where the gas is compressed and the velocity is close 





SiO 


CO(/>14) a 


H 2 0-0^ 




Tkin 








Tkin 




K 


cm -3 


K 


cm -3 


K 


L1448 


Bl 


>500 


8xl0 5 


1200 


6xl0 4 


950-1900 


Rl 


>500 


10 6 


1200 


6xl0 4 


950-1900 


R4 


200 


2.5xl0 5 


600-650 


1-3X10 4 


600-650 


L1157 


Bl 


150-300 


3xl0 5 


350-800 


5-60X10 5 


560-1040 


B2 


200-300 


2xl0 5 


350-800 


5-60X10 5 


560-1040 


R 


50-100 


1-5X10 5 






570-860 



a Derived from ISO-LWS data: Nisini et al. (2000) for L1448 and 
Giannini et al. (2001) for LI 157. 

b Derived from ISO-SWS and ISO-CAM data: Nisini et al. (2000) for 
L1448 and Cabrit et al. (1998) for LI 157. 



to the shock velocity, and lower temperatures occur further 
out where the gas has been slowed down. In L1157-B1, how- 
ever, we find a puzzling result, since most of the gas emission 
comes from the lower excitation component close to the ambi- 
ent velocity, a feature that is not predicted by shock models, as 
we will further discuss in Sect. 4. Jimenez-Serra et al. (2004) 
have suggested that SiO emission at ambient velocity associ- 
ated with outflows could be a signature of the interaction be- 
tween the magnetic shock precursor in a C-type shock, with 
the ambient pre- shock gas. This interpretation in the case of 
LI 157-B1 is however in contrast with the large SiO abundance 
of the low velocity component. In fact, the derived velocity av- 
eraged column density of 8xl0 13 cm~ 2 is almost entirely dom- 
inated by the gas component at ambient velocity. Such a high 
column density implies that this component should have expe- 
rienced significant shock reprocessing. In addition, the ambient 
component due to the precursor is expected to have very narrow 
widths, of the order of 1km s" 1 or less, as observed by Jimenez- 
Serra et al. (2004) in L1448, while in LI 157 the width of the 
component at ambient velocity is ~7km s" 1 . The inferred differ- 
ent excitation conditions of the low and high- velocity compo- 
nents in LI 157 could be instead due to an evolutionary effect, 
where the SiO gas at ambient velocity in LI 157 arises from de- 
celerated post-shocked material, that has already slowed down, 
as suggested in Codella et al. (1999). 

The velocity structure of the LI 448 Bl and Rl bullets, cou- 
pled with the inferred variations of physical parameters with 
the velocity, resemble very closely those found in optical/IR 
jet-knots near the protostar. Indeed, the base of the optical jets 
structure is characterised by the presence of both low and high 
velocity components, with this latter being the brightest and 
having the highest excitation conditions (e.g. Ray et al. 2006). 
This similarity, and the evidence that the collimated EHV bul- 
lets are hotter and denser than usually found in the CO en- 
trained outflow, suggest that they may indeed trace the main jet 
from the embedded protostar. Such an interpretation has been 
also suggested for the EHV SiO jet of HH211 on the basis of 
its excitation and velocity structure derived from SMA obser- 
vations (Hirano et al. 2006). 
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Fig. 5. SiO brightness temperatures predicted by LVG model 
calculations for different temperatures and for densities of 10 5 
cm -3 , (upper panel) and 10 6 cm" 3 (lower panel). A ratio be- 
tween column density and velocity spread Ns/o/AV equal to 
10 13 cm~ 2 km -1 s has been considered in this plot. 



4. Comparison with other warm gas tracers 

In the physical conditions traced by the SiO emission, i.e. 
T -100-500 K and n ~10 5 -10 6 cm" 3 , the main gas coolants 
are pure-rotational H2 lines, CO transitions with J up between 
10 and 20, and H2O rotational lines. H2O cooling is only signif- 
icant if the conditions are appropriate to have large water abun- 
dances; this is indeed the case in low- velocity C-type shocks 
where most of the gas-phase SiO is expected to be produced. 
Most of the cooling from these molecules occurs at mid and 
far IR wavelengths. Since both the investigated outflows have 
been observed with ISO we can compare our diagnostic on the 
SiO emission with the results obtained by ISO on the other 
main coolants (see Table 4). In L1448, strong CO with / >14 
and H2O emission is associated with the Bl and Rl bullets and 
the inferred physical conditions are r=1200 K and ft# 2 =6xl0 4 
cm -3 (Nisini et al. 1999, 2000). The temperature is also con- 
sistent with the H2 0-0 lines detected by SWS. The SiO emis- 
sion seems to be confined to a denser region with respect to the 
other coolant, which implies a large density stratification in- 
side the bullets. In the R4 bullet the conditions derived from the 
ISO-observed CO lines are n Hl ~l-3xl0 4 cm" 3 and T -600 K. 
These are consistent with the lower excitation found in this po- 
sition with SiO, although for this latter a higher density and 
lower temperature are estimated. From the comparison of our 
derived SiO column densities ad the N(CO) estimated from the 
ISO data (~10 17 cm" 2 , Nisini et al. 2000), we obtain a SiO 
abundance of ~10" 7 and 3xl0" 8 in the Bl/Rl and R4 posi- 
tions, respectively (Table 5). 

In LI 157 the temperature and density estimated from the 
ISO data ranges between 350-800 K and 5xl0 5 -6xl0 6 cm" 3 
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Fig. 6. SiO line intensities observed in knots Bl, Rl and R4 
of the LI 44 8 outflow, normalized to the J =5-4 intensity and 
compared with the LVG model computations. Solid lines refer 
to the best fit through the data, while the dashed line represents 
the model fit at the lowest temperature which is still compatible 
with the data errors. The plotted line ratios uncertainty consider 
both the individual rms errors and an absolute calibration un- 
certainty of 10% for the lines with / <8 and 20% for the /=10 
and 1 1 lines 



respectively (Giannini et al. 2001), thus fairly consistent with 
the parameters derived from SiO. Temperatures in the same 
range are also found from the H2 0-0 lines observed with 
ISOCAM (Cabrit et al. 1998). Therefore it appears that the 
physical conditions in the B 1 bullets are more uniform, with- 
out strong gradients. The SiO abundance is, in both Bl and R 
bullets, of the order of 6-8xl0" 8 . 

From the comparison of the SiO column density and the 
column density of the warm water traced by ISO, we derive a 
SiO/H 2 abundance ratio of ~ 2-3xl0" 4 in L1448 and ~10" 3 
in LI 157 (see Table 5). The difference in these two values 
cannot be explained by the different excitation conditions of 
the two flows. For low velocity shocks, the SiO/F^O ratio 
is sensitive to the shock speed, since shocks with velocity < 
20km s" 1 do not efficiently release silicon by sputtering of dust 
grains, while still being very efficient to produce of large H2O 
abundances by means of high temperature (T^300 K) chemical 
reactions. However, both the maximum velocity and the excita- 
tion conditions observed in LI 157 indicate that here the shock 
speed should be lower than in LI 44 8. The different SiO/F^O 
ratio observed could be due instead to an evolutionary effect, 
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Fig. 7. The same as Fig. H but for the knots R, Bl and B2 of 
LI 157. Filled circles refer to observed ratios assuming beam 
filling equal 1 , while open circles are the observed ratios cor- 
rected assuming an emitting size of ~20"for R and Bl, and 
30"for B2 (see Section 3). The best fits obtained considering 
only the / <5 lines (solid) and all the beam filling corrected 
lines (dashed) are reported. 

Table 5. SiO and H2O abundances 
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a Assuming the H 2 abundance derived by Nisini et al. (2000)(L1448) 
and Giannini et al. (2001)(L1157). 



if the two molecules have different destruction timescales af- 
ter the passage of the shock. The timescales for H2O depletion 
onto the dust grains range between 10 4 -10 5 yr depending on 
the density (Bergin et al. 1998). Similar timescales are also 
predicted for the SiO destruction either by depletion on dust 
grains or by reactions with OH which bring to the formation of 
Si0 2 (Codella et al. 1999, Schilke et al. 1998). In conclusion, 
no simple considerations based on excitation conditions and 
evolutionary timescales seem to explain the different observed 
abundance ratio. More detailed chemical models of the two re- 
gions are required for a proper understanding of this point. 




V LSR (km s ) 



Fig. 8. Comparison between the intensity profiles vs velocity 
of the SiO 5-4 and 8-7 lines in the L1448 Rl and Bl bullets 
(upper) and the LI 157 Bl bullet (lower). The dot-dashed and 
grey lines are the individual profiles of the 5-4 and 8-7 lines, 
respectively, while the black line represents the profile of the 
line ratio. There is a trend, more evident for LI 157 Bl and 
L1448 Bl, showing that this ratio increases with the velocity. 



5. Comparison with predictions from shock 
models 

The multi-line analysis performed can be used to check if the 
derived excitation conditions and abundances are in agreement 
with the expectations from shock models. 

The SiO molecule is produced in C-type shocks by the in- 
jection into the gas-phase of silicon by sputtering and/or grain- 
grain collisions, followed by gas-phase reactions with O and 
2 (Schilke et al. [T997l Caselli et al. [T997l . To be efficient in 
producing the large observed column densities of SiO by sput- 
tering, such shocks need to have speeds of at least 25 kms" 1 . 
At lower velocities (Y s ~ 20km s" 1 ) and high density, grain- 
grain collisions may be more efficient than sputtering for the 
production of gas-phase elemental silicon. At such speeds, the 
neutral temperature of the post- shocked gas rises above 1000 
K and then gently declines to return to its pre- shock values in a 
region of ~5xl0 16 cm (~10"at a distance of 350 pc, i.e. com- 
parable with the bullets dimensions). 

We can compare the observed line intensity distributions as 
a function of / (Figs. 5,6), with those predicted by the Schilke 
et al. (1987) shock models as a function of pre-shock density 
and velocity (see their Fig. 6). The line ratios derived for the 
LI 448 bullets are consistent with n(H2) ~ 10 6 cm -3 and shock 
velocities larger than 30 kms" 1 . LI 157 Bl is instead better fit 
with a shock at lower density, between 10 5 and 10 6 cm" 3 . The 
models are not very sensitive to the shock velocity for v^ larger 
than 30 km s" 1 . This is a consequence of the fact that the post- 
shocked temperature is a function of the velocity, and that the 
SiO line ratios become insensitive to temperature variations for 
T larger than -500 K. However, the peak post-shock tempera- 
ture reaches values larger than 1000 K already in a 25 kms" 1 
shock (e.g. Kaufman & Neufeld 1996), thus the SiO emission 
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region is indeed expected to be very warm. Observations of 
submillimetre tracers with higher excitation energies, such as 
the high- J CO lines, are needed to better constrain the temper- 
ature values inside the bullets. 

The fractional SiO abundances derived in Table 5 are also 
consistent with those predicted by models in the regions of the 
post- shocked gas having the larger temperatures (between 10" 8 
and 10" 7 ). In further out regions, where the temperature drops 
close to the pre-shock values, SiO is in part converted to Si02 
by reactions with OH, and its abundance decreases to values 
below 10" 9 . The averaged high observed abundances therefore 
testify that we are indeed observing the warmer regions just 
behind the shock front. 

We can finally compare the observed line profiles with that 
predicted by models. Schilke et al. (1997) computed how the 
profiles change as a function of the rotational quantum number, 
both in the case of sputtering of Si-bearing material in grain 
cores and for the case of mantle sputtering. The overall pro- 
files produced in the two cases are similar to those observed 
in the LI 44 8 bullets: very asymmetric, with a gradual decrease 
of intensity for velocities lower than peak velocity and a very 
sharp decrease on the other side. In these bullets the inten- 
sity peaks occur at a Vlsr velocity of about 60km s" 1 which, 
deconvolved for inclination effects, corresponds to a total ve- 
locity larger than 150 kms" 1 (for an inclination angle of 69°, 
Girard & Acord 2001 ). This is much larger than the estimated 
shock velocity, implying that probably the bullet is moving in 
a medium that has been already put into motion, and thus that 
the shock occurs at the much lower velocity given by the dif- 
ferential speed among the two fluids. 

In contrast with L1448, the LI 157-Bl profiles do not ap- 
pears consistent with the predicted ones, since the brightness 
peak is close to the ambient velocity and the intensity declines 
towards the higher velocity. We have already pointed out and 
discussed this behaviour in the previous sections, suggesting 
that the peak at ambient velocity could be due to a shock that 
has already slowed down but still retaining the original SiO 
abundance. Models taking into accounts time evolution effects 
are needed to test this hypothesis. In addition, geometrical ef- 
fects can be also important and are expected to play a role in the 
line profile appearance. Bow shock models, for example, pro- 
duce line profiles that are highly dependent upon the bow incli- 
nation angle and for some configurations produce shapes simi- 
lar to those observed in the LI 157 SiO line (see e.g. Schultz et 
al. 2005). 

Another effect predicted by the Schilke et al. (1997) model 
is that the line profile broadens for lines with increasing /. Such 
a trend is not seen in our observed lines. In Fig. 8, for example, 
the profiles of the 2-1 and 8-7 lines are compared: in both the 
objects (LI 44 8 and LI 157) the two lines are equally broaden. 

In conclusion, current shock models fit reasonably well the 
average physical conditions and abundances of the SiO gas de- 
rived from our multiline analysis in both the EHV bullets of 
LI 448 and in the lower velocity clumps of LI 157. However, the 
line profiles are still not well reproduced, especially when pro- 
file variations among lines with different excitation are taken 
into account. Effects such as a time-dependent chemistry and 
a geometry different from plane-parallel need to be taken into 



account for a more quantitative comparison with the observed 
profiles. In addition, Interferometric observations should be 
able to resolve spatially the emitting regions contributing to the 
different components, thus allowing us to constrain the geom- 
etry to be adopted in a more detailed modelling. It would be 
also instructive to compare the observed SiO profiles with the 
profiles of other molecules similarily affected by a time depen- 
dent chemistry. H2O, in particular, will be studied in details 
by means of the Herschel observatory, providing valuable con- 
straints on the shock model chemistry. 
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Fig. 9. Comparison between the intensity profiles of the SiO 2- 
1 (black) and 8-7 (gray) lines in the L1448 Rl and Bl bullets 
(upper) and the LI 157 Bl bullet (lower). The L1448 8-7 line 
profile toward the mm source central position has been con- 
volved to the beam of the 2-1 transition in order to compare 
the same emitting area. 



6. Conclusions 

In this paper we present the first SiO multi-line analysis (from 
7=2-1 to J=l 1-10) of the molecular bullets along the outflows 
of the Class sources L1448-mm and L1157-mm. In L1448- 
mm, the bullets represent collimated and regularly spaced 
clumps travelling at velocities of hundreds of kms" 1 . In con- 
trast, the LI 157 bullets have lower velocities and are charac- 
terised by a rich chemistry. From our analysis, we investigated 
the physical properties of these different types of bullets, ex- 
ploring their excitation conditions and the capability of current 
models to reproduce the various SiO line ratios and profiles. 
The main results of our analysis are the following: 

- In L1448, the SiO line ratios of the EHV molecular bul- 
lets Bl and Rl close to the driving source are consistent 
with a gas at n Hl ~10 6 cm" 3 and ^ 500 K. Conversely, 
the bullet R4, located at the apex of the red lobe molecular 
outflow and very strong in the SiO 2-1 emission, is asso- 
ciated with lower excitation gas (n Hl - 2.5 10 5 cm" 3 and 
Tkin = 200 K). Also the low velocity bullets in LI 157 show 
lower excitation conditions, with hh 2 -1-5x10 cm J and 
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Tkin ~ 100-300 K. In both the sources, there is however a 
clear evidence of the presence of velocity components hav- 
ing different excitation conditions, with the denser and/or 
warmer gas associated with the gas at the higher speed. 

- While in LI 448 the bullets emission is dominated by the 
highly excited and high velocity gas, in LI 157 the peak of 
the emission at the ambient velocity is associated with the 
lower excitation conditions. Such a difference could be due 
to evolutionary effects, where the LI 157 SiO gas at ambient 
velocity is due to an older shock that has already slowed 
down and is progressively reaching the physical conditions 
of the ambient medium. 

- The physical conditions averaged over all the velocity com- 
ponents are consistent with those inferred by means of mid- 
and far-IR lines observed with the ISO satellite. In LI 448, 
however, the high- J CO lines detected by ISO-LWS trace a 
gas with a density about an order of magnitude lower than 
that measured with the SiO lines, indicating that gradients 
in the physical parameters can be present in the bullets re- 
gion. 

- The column density of the SiO gas does not show large vari- 
ations in the different bullets, ranging between 2xl0 13 and 
10 14 cm" 2 . A comparison with the ISO-LWS derived abun- 
dance of H2O , a molecule also very sensible to the shock 
conditions, shows that the SiO/F^O abundance ratio is ~ 
2-3xl0" 4 in L1448 and ~1(T 3 in LI 157. Such a different 
ratio is not easily explained by either the local excitation 
conditions and the different flow evolutionary timescales. 

- Finally, we find that current shock models for the SiO ex- 
citation reproduce quite well the observed line ratios when 
shocks with velocities larger than ~ 25km s" 1 are consid- 
ered. The models however fail to predict the observed pro- 
file variations as a function of the line excitation, as well as 
the strong emission at ambient velocity observed in LI 157. 
Effects such as a time-dependent chemistry and a geometry 
different from plane-parallel need to be taken into account 
for a more quantitative comparison with the observed pro- 
files. 

- We also point out that the velocity structure and the vari- 
ations of physical parameters with the velocity observed 
in the EHV bullets in LI 448 resemble very closely those 
found in optical/IR jet-knots near the protostar, suggesting 
that the mm collimated bullets have a similar origin and 
excitation mechanism. 
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